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ABSTRACT 

We illustrate the possibility of light trapping and tunneling in periodic arrays of metallic nanopar- 
ticles. A controllable minimum in the transmission spectra of such constructs arises from a collective 
plasmon resonance phenomenon, where an incident plane wave sharply localizes in the vertical direction, 
remaining delocalized in the direction parallel to the crystal plane. Using hybrid arrays of different struc- 
tures or different materials, we apply the trapping effect to structure the eigen-mode spectrum, introduce 
overlapping resonances, and hence direct the light in space in a wavelength-sensitive fashion. 



; 1 Introduction 

■ Tremendous progress in fabrication methods of metallodielectric nanostructures,[l] on the one 
. hand, and rapid development of laser technology, [2j on the other, have recently enabled a wide 

variety of exciting opportunities for light manipulation in subdiffraction length scales. These range 
\ from understanding of fundamental concepts of plasmon-polariton dynamics, [3] through coherent 

and optimal control of electromagnetic energy propagation via metallic structures, [1] to different 

applications in modern optical device technology, [5] single atom-molecule manipulations, [6] and 
. ultrasensitive detection of biological molecules. [7] 

The physical basis of these applications is the surface plasmon resonance phenomenon, owing 

to which nobel metal particles of nanosize strongly enhance electromagnetic (EM) fields tuned 

■ to resonance with the frequency of collective oscillation of conductive electrons confined in the 

■ nanoparticle. It follows from the Mie theory, [8] (and was seen numerically as well as experi- 
. mentally) that the resonance wavelength depends crucially on the particle size and shape. This 
'. sensitivity is the origin of many of the applications of nobel metal nanoparticles, including their 
\ use as sensors, [9l HO] medical diagnostics, and standards. [12] 

Like the response to incident light of single nanoparticles (NPs) , the collective response of NP 

■ chains has been the topic of significant experimental [13] and theoretical interest in recent years. [Ti] 
Particularly interesting is the possibility of guiding light in the nanoscale via metallic NP chains 

, and junctions, Here, one excites the construct locally in space at one end of the chain, typically 
[ by means of a nanosized tip, and detects the transmitted energy, again with spatial resolution, at 
\ the other end, e.g., by means of a dye molecule. Much progress on this question not withstanding, 

■ the adverse effect of losses have been identified as the major hurdle to realizing waveguides based on 
. NP chains. [3] Films[16j and wires[T7| suffer much less from losses, but lack the local enhancement 
! and mode specificity of NPs. 

Similarly fascinating, are the properties and applications of plasmonic nanocrystals. By nanocrys- 
tals, we refer to ordered two- or three-dimensional arrays of nanoparticles, periodic in the lateral 
direction, that are collectively excited by an incident (delocalized) plane wave. Metallodielectric 
crystals exhibit omnidirectional band gaps[T8] along with a broadband absorption. [19j Experi- 
mental measurements of optical properties of one- [20] and two-dimensional [2 Ij periodic metal 
nanostructures have been successfully performed. We also note that, not only periodic struc- 
tures composed of NPs exhibit interesting phenomena but also periodic arrays of sub-wavelength 
holes have surprising optical transmission. [22] The physics underlying such constructs has aspects 
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in common with that underlying finite NP chains and junctions, it rehes on the collective re- 
sponse of ordered NPs that are sufficiently closely-spaced to interact via dipole and quadruple 
forces. The new features arise from the periodicity in one direction, which gives rise to certain 
solid state features, such as a band structure, and justifies the name "crystals". The practical 
interest in nanoplasmonic crystals steams from the premise of scaling to the nano-domain sev- 
eral of the attractive properties of micron-sized analogs, such as photonic crystals and band-gap 
materials, which have been extensively studied experimentally and theoretically for their many 
applications. [23] Of particular interest are properties such as hght-induced transparency. [2^ self- 
collimation,[10j and transmission through sharp bends. [25] While the essential physics responsible 
for these features in the macro- world is quite different from that underlying the nanoscale variant, 
the effects, as illustrated below, are similar. 

We show in the following sections that plasmonic nanocrystals offer also new and fundamentally 
interesting phenomena, that are unique to the nanosize and do not have a direct analog in their 
micron-sized counterparts. In particular, we illustrate the possibilities of introducing overlapping 
resonances in the eigen-mode spectrum as a control tool, and of funneling the incident plane wave 
in a predetermined direction with wavelength selectivity. In the next section we brieffy outline the 
theory and the numerical approach taken. Section [3] presents and discusses our results and the 
final section concludes, with an outlook to future research. 



2 Theory 

The interaction between the light and the metal nanoconstructs is simulated using a finite- 
difference time-domain approach. p6] We restrict attention to two-dimensional structures, and 
consider the TE^ mode (transverse-electric mode with respect to z) of the EM field. Maxwell's 
equations in this case read, 
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where E^., Ey, and Hz denote the Cartesian components of the EM field, Jx and Jy are the Cartesian 
components of the current density, J, Eejf is the effective dielectric constant defined below, and fiQ 
is the magnetic permeability of free space. The metal structures are described within the Drude 
model with a complex valued, frequency-dependent dielectric constant. 



e (uj) = Uoo o ^^.-p ) > (2) 



where Eq is the electric permittivity of free space, £00 is the dimensionless infinite frequency limit 
of the dielectric constant, Up is the bulk plasmon frequency, and F is the damping rate. In the 
numerical simulations we use the following two sets of parameters, which describe fairly well silver 
and gold nanoparticles in the incident wavelength regime of interest. For silver nanoconstructs, 
Eoo = 8.926, LUp = 1.7601 x 10^^ rad/sec, and F = 3.0841 x 10^^ rad/sec. For gold nanoconstructs, 
£00 = 9.84, Up = 1.3673 x 10^^ rad/sec, and F = 1.0179 x 10^^ rad/sec. 

In regions of space occupied by metallic particles, the material dispersion gives rise to time- 
dependent current and Maxwell's equations ([1]) are supplemented by the additional set of equations 
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for the current density J.[28| [29] 

(97 

— = aJ + hE, (3) 

where a = —T, b = e^Up and the effective dielectric constant in Eq. ([T]) is Seff = ^o^oo • In 
the surrounding free space Eeff = Eq and a = b = 0. In order to avoid nonphysical reflections of 
outgoing electromagnetic waves from the grid boundaries, we implement perfectly matched layers 
(PML) of absorbing boundaries [30j of depth of 16 spatial steps each. 

As an initial condition we employ a plane wave generated along a line, each point of which is 
driven by a time-dependent function of the form, 

Ei^c{t) = e^E{t) cos ut, (4) 

where is a unit vector along the x-axis, E (t) = i^osin^ {.T^t/r) is the pulse envelope, r is the 
pulse duration {E {t > r) = 0), and uj is the optical frequency. In all simulations r = 50 fs, the 
total propagation time is 100 fs, with a spatial step size of Sx = 6y = 1.5 nm and a temporal step 
size of 6t = 6x/ (1.5c), where c denotes the speed of light in vacuum. All simulations have been 
performed on distributed memory parallel computers at the National Energy Research Scientific 
Computing Center and San Diego Supercomputer Center. The parallel technique used in our 
simulations is described in detail in Ref. 31. 



3 Results and Discussion 

We open this section by exploring the properties of the simplest plasmonic crystal that could be 
envisioned. This construct will serve to introduce the basic phenomena that underlie several of the 
much more complicated constructs discussed below. The nano-crystal envisioned is schematically 
depicted in the inset of Fig. 1. It is composed of a periodic array of silver NPs with 55 nm 
diameter and 66 nm center-to-center distance. The NPs are excited by an x-polarized plane wave 
propagating in the negative ^/-direction, and the EM energy is detected along the horizontal line 
shown in dashed in the inset. The transmission through the nanostructures is expressed in terms 
of the ratio of the time-averaged EM energy in the presence of the periodic array, W^totah to that 
in its absence, Wi^c- This ratio is shown in the main frame of Fig. 1 for five arrays with different 
number of layers as a function of the incident wavelength. The transmission is close to perfect 
throughout the wavelength regime of relevance, safe for a narrow window centered at A ~ 360 nm, 
where it drops to a sharp minimum. The transmission minimum is sharper the larger the number 
of layers, converging to essentially zero transmission (VFtotai/W^inc ~ 2 x 10"'') with 5 layers. 

The origin of the transmission minimum is the excitation of longitudinal (x-polarized) plasmons, 
due to which an incident light propagating along the y-axis bends as it enters the periodic array 
and is guided along the horizontal axis. The efficiency with which the light takes the sharp corner, 
and the spatial localization of its subsequent path are investigated below. Our results provide 
a likely explanation to the finding of Ref. 19, where a strong broadband absorption in metallic 
photonic crystals is reported. As shown below, the resonance feature is due to an interesting 
trapping phenomenon, whereby the electromagnetic energy is localized along the surface of the 
crystal. Such collective waves that are fully delocalized laterally but sharply localized in the 
vertical direction are familiar in electronic systems, where they arise from image potentials. [32] 
We are not aware of previous studies that have observed such phenomena in plasmonic waves. 

Both the broad transmission window and the ability of the crystal to re-direct the propagating 
EM energy by 90 degrees and guide it in parallel to the crystal plane at the transmission mini- 
mum wavelength are of significant practical interest, as they address two long researched goals of 
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nanophotonics. [3] More interestingly from a fundamental perspective, the resonance phenomenon 
that underlies the effect offers exciting opportunities for optical control. Such resonance features 
have been intensively studied in atomic and molecular physics, where their manifestation in radia- 
tive lineshapes is conveniently described within the Feschbach partitioning of space into a bound 
and a scattering manifold. The present system can be analyzed using the machinery of Fano and 
Feschbach by analogy to the well-studied case of matter waves, but offers the advantage over the 
latter case that the interaction coupling the discrete state to the continuum can be systematically 
tuned by design of the array. 

We proceed by exploring the dependence of the resonance feature on the experimentally variable 
parameters of the array. In Fig. 2 we focus on the single layer array and gradually modify the 
resonance properties by coating the silver nano-spheres by a non-dispersive insulator. The inset 
of Fig. 2 depicts schematically the array envisioned. We remark that core-shell nanoparticles of 
the sort shown in the inset of Fig. 2 have been synthesized in several recent studies that have 
noted the potential applications and interesting properties of such constructs. Several trends in 
Fig. 2, and in related studies that are not shown here, can be qualitatively understood from the 
Mie expression for the extinction cross section of a single sphere. Both an increase in the dielectric 
constant of the core material and an increase in the dielectric constant of the shell material red- 
shift the resonance position. An increase in the shell dielectric constant has the additional effect 
of broadening the resonance, see Fig. 2. 

Of particular interest are devices that support multiple resonances. In molecular systems, the 
case of overlapping resonances has been studied extensively for its fundamental interest as well as 
for intriguing opportunities for coherent control. [33] In the nanoplasmonics considered here, both 
the case of isolated and that of overlapping multiple resonances are of fundamental and practical 
interest, and both can be readily realized by proper design of the nano-crystal. Figure 3 explores 
both the case of non-overlapping and that of overlapping resonances. To that end we consider 
hybrid nanocrystals consisting of ordered layers of two types of metallic nanoparticles and include 
both spheres and ellipsoids. The insets depict the two nano-crystals considered, one of which (inset 
a) combines silver and gold nano-spheres and the other (inset b) considers alternating layers of 
spheres and ellipsoids, the top two of which are silver, and the bottom two gold. The combination 
of Au and Ag NPs provides a model for study of non-overlapping resonances, since the plasmon 
resonance wavelengths of Ag and Au are separated by ca 120 nm. The combination of spheres and 
ellipsoids provides a model for study of overlapping resonances, since the shape-dependence of the 
resonance wavelength is small, leading to a shift of the line center comparable to the resonance 
width. 

The main frame of Fig. 3 shows simulations of the ratio, H^totai/W^inc, for the two hybrid 
constructs depicted in the insets, with the solid curve corresponding to the bimetallic nano-crystal 
of inset (a) and the dashed curve corresponding to the mixed crystal of inset (b). Both spectra 
show the expected well-defined double-peak structure. The spectrum generated by construct (b) 
exhibits also significant broadening and additional structure, arising from the partially overlapping 
resonances of the spheres and ellipsoids. 

The results of Fig. 3 can be understood through an eigenmode analysis of the elements con- 
stituting the nanocrystals. The eigenmodes of a given construct are computed by excitation with 
a broad band, short (0.36 fs) pulse, followed by long time (2 ps) simulation of the x-component 
of the resultant scattered electric field. Fourier transform of the field yields the power spectrum 
/(A). The solid and dashed curves in Fig. 4 correspond to the cases of single (periodic) layers 
of silver and gold nano-spheres, respectively, whereas the dot-dashed curve corresponds to one 
Ag and one Au nano-spheres layer and the dotted curve to one Ag spheres and one Ag ellipsoids 
layer. Several interesting features are noted, including Fano-type interference between direct and 
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resonant transmission resulting in marked asymmetry (dashed curve) , non-overlapping resonances 
(dot-dashed curve) and overlapping resonances (dotted curve). 

The plasmon resonances of silver and gold linear periodic chains are separated by almost 120 
nm, as illustrated in Fig. 4, resulting not only in selective transmission of light by the hybrid 
nano-structure but also in separation of different resonant wavelengths in space. For example, 
wavelengths near the plasmon resonance of silver chains, A ~ 380 nm, are trapped in the first two 
silver layers (see the insets in Fig. 3), whereas wavelengths near the plasmon resonance of gold 
chains (A ~ 500 nm) are transmitted through the silver layers and excite the bottom two gold 
chains. The wavelength-dependent transparency of different parts of a nano-structure is due to 
the fact that the effective thickness of the skin-layer at the wavelengths considered is of the same 
order of magnitude as the size of the NPs. The skin-depth {d) at optical frequencies (u) much 
larger than the plasmon frequency {oOp) can be estimated[3lj as ~ c/ujp, which leads to 17 nm 
and 22 nm for silver and gold, respectively. 

The wavelength-sensitive transparency of different parts of nanostructures suggests exciting 
opportunities for use of non-symmetric (with respect to y) constructs for light manipulation and 
directed transport. Two examples of light guidance are presented in Fig. 5, where we simulate 
hybrid structures composed of silver and gold NPs that are finite in both x- and y-directions. At 
A = 340 nm, for instance, the left part of the nano-crystal depicted in the inset (a) of Fig. 5 traps 
the light, while the right part, composed of gold NPs is transparent. Excited at their plasmon 
resonance wavelength, silver layers guide the EM energy to the right. The structure shown in inset 
(b) of Fig. 5, supports both horizontal directions for EM energy guiding. This is illustrated in Fig. 
6, where we plot the spatial time-averaged distribution of the ratio, Wtotai/W^inc, at the incident 
wavelength corresponding to the plasmon resonance of silver NPs. Clearly, an incident plane wave, 
initially symmetric, is broken by the hybrid structure into two parts, one of which is guided to the 
left, and the other to the right. 

In summary, we illustrated selective trapping and guidance of light by periodic metallic nanocrys- 
tals. Owing to the excitation of collective longitudinal plasmons in linear chains of closely spaced 
NPs, an incident wave can be selectively localized in the direction vertical to the light propaga- 
tion, remaining delocalized laterally. By proper device design, EM excitation can be thus efficiently 
guided in space. A variety of extensions of the phenomena of selective light trapping and funneling 
can be envisioned. Three dimensional nanocrystals could be used to guide EM energy in 3D space. 
The interaction of white light with periodic arrays could serve to make interesting optical sources. 
Likewise inviting is the possibility of guiding light in circular or helical pathways using properly 
structured plasmonic nanoconstructs. 
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Figure 1: The ratio of time- averaged EM energy calculated in the presence and in the absence of 
the nanoconstruct schematically depicted in the inset as a function of the incident wavelength, A. 
The solid curve shows the ratio for a single horizontal layer of NPs, the dashed curve for two layers, 
the dotted curve for three layers, the dash-dotted curve four layers, and the dash-dot-dotted curve 
for five layers. The inset shows the geometry envisioned, where the dashed vertical lines represent 
periodic boundaries in the lateral direction. The nanoconstruct is excited by an x-polarized plane 
wave propagating in the negative y direction and the EM energy is calculated along the horizontal 
line shown as a dashed line beneath the structure in the far-field zone (the distance of the detection 
line from the crystal is 600 nm). 
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Figure 2: As in Fig. 1, for a single layer of silver NPs coated by an insulator (see inset). The 
solid curve shows the ratio PVtotai/W^inc for bare particles (no coating, corresponding to e = 1), 
the dashed curve for coating dielectric constant of e = 1.5, the dotted curve for e = 2.5, and the 
dash-dotted curve for e = 3.5. The silver NP diameter is 35 nm, the diameter of the dielectric 
coating layer is 55 nm, and the center-to-center distance is 66 nm. 
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Figure 3: The ratio Wtotai/W^inc as a function of the incident wavelength for the two nanoconstructs 
shown schematically in the insets, (a) A silver-gold hybrid structure composed of spheres of 
diameter 32 nm and center-to-center distance 42 nm. (b) A silver-gold hybrid structure composed 
of spheres and ellipsoids with major (minor) axes of 20 (32) nm and center-to-center distance of 
42 nm. Gold NPs are shown as black, silver particles as grey spheres. The solid curve in the main 
frame corresponds to structure (a) and the dashed curve to structure (b). 
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Figure 4: Power spectrum, |/ (A)|, of the x-component, E^, of the scattered electric field in Fig. 3 
as a function of the incident wavelength, A. The solid curve corresponds to a single layer of silver 
NPs, the dashed curve to a layer of gold NPs, the dash-dotted curve to interacting layers of silver 
and gold NPs, and the dotted curve to silver layers of spheres and ellipsoids. 
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Figure 5: The ratio Wtotai/W^inc as a function of the incident wavelength for the two finite nanocon- 
structs shown schematically in the insets. Silver NPs are shown as grey spheres and gold NPs as 
black spheres. The solid curve in the main frame presents simulations for structure (a), and the 
dashed curve for structure (b). Both silver and gold particles are 40 nm in diameter and the 
center-to-center distance is 50 nm. 
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Figure 6: Time-averaged ratio, VFtotai/W^inc, as a function of x and y for the structure shown in 
inset (b) of Fig. 5 at an incident wavelength of A = 340 nm, corresponding to the excitation of 
longitudinal plasmons in silver NPs. 
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